arXiv:1501.02335v2 [quant-ph] 22 Mar 2015 


Characterizing non-Markovianity via quantum interferometric power 


Himadri Shekhar Dhar , 1 - 2 Manabendra Nath Bera , 1 and Gerardo Adesso 3 

1 Harish-Chandra Research Institute, Chhatnag Road, Jhunsi, Allahabad 211 019, India 

2 School of Physical Sciences, Jawaharlal Nehru University, New Delhi 110 067, India 
3 School of Mathematical Sciences, 

The University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom 

Non-Markovian evolution in open quantum systems is often characterized in terms of the backflow of in¬ 
formation from environment to system and is thus an important facet in investigating the performance and 
robustness of quantum information protocols. In this work, we explore non-Markovianity through the break¬ 
down of monotonicity of a metrological figure of merit, called the quantum interferometric power, which is 
based on the minimal quantum Fisher information obtained by local unitary evolution of one part of the sys¬ 
tem, and can be interpreted as a quantifier of quantum correlations beyond entanglement. We investigate our 
proposed non-Markovianity indicator in two relevant examples. First, we consider the action of a single-party 
dephasing channel on a maximally entangled two-qubit state, by applying the Jamiolkowski-Choi isomorphism. 
We observe that the proposed measure is consistent with established non-Markovianity quantifiers defined us¬ 
ing other approaches based on dynamical divisibility, distinguishability, and breakdown of monotonicity for the 
quantum mutual information. Further, we consider the dynamics of two-qubit Werner states, under the action of 
a local, single-party amplitude damping channel, and observe that the nonmonotonic evolution of the quantum 
interferometric power is more robust than the corresponding one for entanglement in capturing the backflow of 
quantum information associated with the non-Markovian process. Implications for the role of non-Markovianity 
in quantum metrology and possible extensions to continuous variable systems are discussed. 
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I. INTRODUCTION 


The fundamental essence of any information protocol resides 
in the amount of accessible information in a given physical 
system. In a realistic quantum process, the system is rarely 
isolated and commonly interacts with an environment, which 
results in the information being scattered in the typically large 
Hilbert space of the environment. In most theoretical mod¬ 
els, an open system is studied through the dynamics of the 
reduced density matrix, upon tracing over the environmental 
degrees of freedom [1, 2]. The dynamics of the system is 
described by completely positive semigroup maps, or equiva¬ 
lently by the solution of a master equation in Lindblad form 
[3, 4]. This formalism assumes weak system-environment 
coupling, short environment correlation time and a “memory¬ 
less” transfer of information from the system to the environ¬ 
ment leading to information erasure [1,5]. Such a dynamical 
model for open quantum systems is called Markovian. How¬ 
ever, it is readily observed that the Markovian formalism is 
not always optimal or justified when dealing with many im¬ 
portant open quantum systems, especially in complex biolog¬ 
ical models or interacting many-body systems in condensed 
matter physics [6]. More accurately, the system-environment 
interaction needs to be treated as non-Markovian, which tends 
to deviate from the completely positive semigroup dynamics 
[6, 7], thus making the corresponding mathematical formal¬ 
ism difficult. Incidentally, non-Markovian dynamics are not 
“memoryless” and allow a backflow of information from the 
environment [8, 9] (cf. [10]) to the system, a fact which has in¬ 
teresting ramifications from the perspective of quantum infor¬ 
mation theory [11]. For example, non-Markovian processes 
have been shown to preserve entanglement [12] in many-body 
[13] and biomolecular [14] systems, and have been exploited 


in quantum key distribution [15], enhancing precision in quan¬ 
tum metrology [16], and implementing certain quantum in¬ 
formation protocols [17, 18]. Non-Markovianity also plays 
a detrimental role in quantum Darwinism, thus impeding the 
emergence of classical objectivity from a quantum world [19]. 

With recent development of experimental techniques to en¬ 
gineer and control system-environment interactions [20] (see 
also Ref. [1]), there is considerable interest in characterizing 
and quantifying non-Markovian dynamics and investigating 
possible applications in scalable quantum technologies [15- 
18] that are robust against environment-induced decoherence 
[13] or phenomena such as entanglement sudden death [21]. 
Although the concept of non-Markovianity is well established 
in the classical realm [22], its quantum extension is often rid¬ 
dled with inconsistency and subtle variations. This has led to 
a substantial amount of literature attempting to quantitatively 
characterize non-Markovianity based primarily on the non¬ 
monotonic time evolution of some quantum information mea¬ 
sure (for reviews, see [23-25]). Such nonmonotonic behavior 
arises from the nondivisibility of the completely positive and 
trace preserving (CPTP) maps [6] that describe the dynamics 
of the open quantum system, which is perhaps the most es¬ 
tablished marker of non-Markovianity [24-26] (cf. [27]). The 
nondivisibility of a CPTP map is necessary for the occurrence 
of information backflow from the environment or the presence 
of the environment memory [8]. 

A number of non-Markovian measures and witnesses have 
been proposed. Among the most important ones, let us men¬ 
tion those based on the deviation of the dynamical maps from 
divisible CPTP maps [6, 7] and those based on the nonmono¬ 
tonicity of the trace distance or distinguishability [8], entan¬ 
glement [6], quantum mutual information [28], and channel 
capacities [18]. Other significant attempts to quantify non- 
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FIG. 1. (Color online). A scheme showing the protocol under con¬ 
sideration. The system a interacts with the environment, while the 
quantum interferometric power is computed by applying local uni- 
taries on the ancilla part b, which serves as the measuring apparatus. 


QIP is qualitatively consistent with measures based on distin- 
guishabilty, divisibilty, and quantum mutual information. In 
Sec. IV, we consider a single-qubit amplitude damping model 
and investigate the flow of QIP in the non-Markovian regime. 
We observe that the measure appropriately captures the back- 
flow of information and is more robust compared to a non- 
Markovianity measure based on entanglement. We discuss 
the results, possible extensions, and potential benefits of the 
introduced non-Markovianity indicator in Sec. V. 

II. CHARACTERIZING NON-MARKOVIANITY VIA 
QUANTUM INTERFEROMETRIC POWER 

A. Quantum interferometric power 


Markovianity include the flow of quantum Fisher informa¬ 
tion [29], fidelity between dynamical time-evolved states [30], 
distinguishability in Gaussian quantum systems in terms of 
fidelity [31], volume of Gaussian states [32], backflow via 
accessible information [33], and local quantum uncertainty 
[34]. Recent proposals have also been made to character¬ 
ize non-Markovianity in direct analogy to entanglement the¬ 
ory [35] and to study the links between system-environment 
quantum correlations and non-Markovian behavior [36]. In¬ 
terestingly, although the different non-Markovian measures 
and witnesses emanate from the dynamical divisibility crite¬ 
ria, the inverse implication is not always true, which makes 
them incompatible with each other for general open system 
dynamics [23, 24, 37] (cf. [28, 38]). 

In this work, we propose to characterize the non- 
Markovianity of an open system evolution through the non¬ 
monotonic behavior of a quantum metrological figure of 
merit, called the quantum interferometric power (QIP) [39, 
40], which is defined in terms of the minimal quantum Fisher 
information [41] obtained by local unitary evolution of one 
part of the system. The QIP is an important information- 
theoretic tool that also quantifies discordlike quantum corre¬ 
lations in a bipartite system [39, 40] and is related to the min¬ 
imum speed of evolution of a quantum system in the projec¬ 
tive Hilbert space [40]. To capture the non-Markovianity in 
open quantum evolutions, we consider a single qubit (say, a) 
as the principal system, interacting with an environment. A 
second qubit (say, b ) plays the role of an ancilla. We con¬ 
sider the action of the environment on the system in terms 
of the dephasing and amplitude damping channels. Using 
the Jamiolkowski-Choi isomorphism [42], single-qubit oper¬ 
ations can be used to study the bipartite (system + ancilla) 
behavior. The QIP of the system is measured by applying 
local unitaries on the ancilla, which acts as a measuring ap¬ 
paratus. The non-Markovianity of the evolution is character¬ 
ized by quantifying the nonmonotonic behavior of the QIP. 
The paper is organized as follows. In Sec. II, we briefly 
present the concept and definition of the QIP and introduce 
a non-Markovianity measure based on its nonmonotonic evo¬ 
lution. In Sec. Ill, we consider a prototypical single-qubit 
dephasing model for the two-qubit (system + ancilla) state, 
and show that the non-Markovianity measure derived using 


The QIP is a metrological figure of merit that quanti¬ 
fies the guaranteed precision enabled by a bipartite probe 
state for the task of black-box quantum parameter estimation 
[39, 40, 43, 44]. Let us consider a bipartite, system + ancilla 
state, p a b, such that the ancilla ( b ) is subject to a local uni¬ 
tary evolution. In this picture, the ancilla acts as a measuring 
device for any operation performed on the system (a) (for an 
illustration, see Fig. 1). The system + ancilla Hamiltonian 
is given by 'H - I„ <8> Hb, where '74 is the local Hamilto¬ 
nian acting on b, and I a is the identity operator acting on a. 
For any bipartite state, p a b and local Hamiltonian 'Hi,, the op¬ 
timal available precision for the estimation of a parameter 0 
encoded in the local unitary Ub - exp(-;0'74X is governed 
by the quantum Fisher information [41], as derived using the 
Cramer-Rao bound [45], which is defined as follows. 

For a quantum state, written in its spectral decomposition 
as Pab = E m e m |0mX0m[, where e m > 0 and £ m e m = 1, the 
quantum Fisher information associated with the local evolu¬ 
tion generated by I a ® Hi, can be written as [41], 




-« z 

m,n: 

e m +en>0 


^ e " ) 2 \{<t>m\Ia ® 'Hb\(pn)\ 2 - 

6 m “ 1 “ &n 


(i) 


The above expression can be equivalently rewritten as follows, 

T(p ab , 'Hb) = ^r(pH 2 b ) - V -^H\((f> m \Ia ® 'Hb\cpn)\ 2 - 

e ™ + e « 

e m +e n > 0 

( 2 ) 

If the generator of the local evolution on the ancilla is not 
known a priori, as in the black-box paradigm for quantum 
metrology [39, 43], then the guaranteed precision enabled by 
the state p a b is given by the QIP (Q), defined as the minimum 
quantum Fisher information over all local Hamiltonians '74 
of a fixed spectral class (a canonical choice is to consider the 
minimization to run over all '74 with nondegenerate, equis- 
paced eigenvalues) [46], namely 

Qipah) = ] -mtr(pab,'H h ), (3) 

4 ‘tti. 


where the | factor is a convenient normalization [39]. 
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In the relevant case where the system a has arbitrary di¬ 
mension, while the ancilla b is a qubit, the choice of local 
Hamiltonians is reduced to < Hb = r ■ &, where \f\ - 1 and 
cr = {cr x , a y , cn] is the vector of the Pauli matrices [39, 40]. 
For such local Hamiltonians, the minimization in Eq. (3) can 
be performed analytically, so that the QIP is computable in 
closed form and given by the expression [39, 40], 

QiPab) = 1 - AZ ax , (4) 

where 4""“ is the highest eigenvalue of the real symmetric 
matrix W with elements [39, 40], 

Wij= V ~~ g '” e '' {(/) m \I a ® er' b \il>„)((/) n \I a ® cr J b \(f> m ). ( 5 ) 

m,n : e '» + e » 

e m +e n >0 

In general, the QIP is a bona fide measure of bipartite quan¬ 
tum correlations beyond entanglement, of the so-called dis¬ 
cord type (see [47] for a review), in the quantum state p a b. 
Namely, Q(p a b) is known to vanish for states with zero dis¬ 
cord from the perspective of subsystem b (known as quantum- 
classical states), is invariant under local unitary operations, 
and reduces to an entanglement monotone for pure quantum 
states. Most importantly for the aims of the present paper, the 
QIP is a monotonically decreasing function under the action of 
arbitrary local CPTP maps on the system a [39, 40]. Further¬ 
more, the QIP can be interpreted as the minimal global speed 
of evolution for the state p c ,b under all local unitary transfor¬ 
mations on the ancilla b, as a consequence of the connection 
between the quantum Fisher information and the Bures metric 
[40, 48, 49], 

The evaluation of the QIP remains computationally 
tractable for higher-dimensional ( d a x d/,) systems, although 
a closed analytical form may not be available. The problem 
can be recast in the form of a minimization of the Hamiltonian 
with respect to a finite number of variables spanning a com¬ 
pact space. This follows by noting that the unitary evolution, 
corresponding to < Hb acting on the ancilla b, can be chosen 
within the special unitary group, without any loss of general¬ 
ity. Furthermore, the QIP can also be reliably computed for 
two-mode Gaussian states in the continuous variable regime 
[43, 44], 


B. Characterizing non-Markovianity 

Let us consider an open quantum system undergoing an 
evolution given by the time-local master equation. 


Here A, -(f) are the Lindblad operators, and y,(f) is the time- 
dependent relaxation rate. 

The quantum evolution is Markovian when y,(f) > 0 for 
each instant of time, V 1 > 0. The dynamical quantum pro¬ 
cess can be then defined in terms of time-ordered CPTP maps, 
such that fl(t 2 , t \) = T exp [ dt' £,(t ')^, where T is the time¬ 
ordering operator. The map Q(t, 0) represents the evolution of 
the system from an initial state (t = 0) to a state at time t. Im¬ 
portantly, such CPTP map satisfies the divisibility criteria, in 
the sense that it can be written as a composition of other time- 
ordered CPTP maps, such that Q(t + dt, 0) = il(t + dt, 0). 

Conversely, for instances where y,(f) < 0, the correspond¬ 
ing dynamical map O (t + dt,t) may not be CPTP and the di¬ 
visibility property of the overall CPTP dynamics is violated. 
The non-divisibility of the dynamical maps given by a time- 
local master equation, of the form in Eq. (6), is the essential 
marker of non-Markovian dynamics [24—26] (cf. [27]). As 
most quantum information quantities are monotonic under lo¬ 
cal CPTP maps, any observation of a nonmonotonic behav¬ 
ior in some reference quantity can be exploited to capture the 
breakdown of Markovianity. 

In this work, we consider the QIP ( Q(p a b )) as our refer¬ 
ence information-theoretic figure of merit, and witness non- 
Markovianity in terms of the nonmonotonicity of the QIP 
under local evolutions of the system a, since for Markovian 
dynamics it must hold that Q(Q. a p a b ) < Q(p a b ), for all local 
CPTP maps Q„ acting on the system a. 

More precisely, to characterize non-Markovianity, let us be¬ 
gin with a basic description of the open system under consid¬ 
eration. Let a be the system interacting locally with an exter¬ 
nal environment and b be the ancilla of the overall bipartite 
state p a b. We calculate the QIP in the bipartite system us¬ 
ing the measure Q(p a b ) defined in Eq. (3), where the ancilla b 
acts as the measuring apparatus on which the local unitaries 
are applied (see Fig. 1). Now we consider the dynamics of 
the system a as described by the time-local master equation 
of Eq. (6). For a Markovian dynamics, the evolution is given 
by a divisible CPTP map Cl a {t). Using the Jamiolkowski-Choi 
isomorphism, the composite dynamics of the overall system, 
Pab(t), is given by 

Pab(t) = (fl a (t) ® Ib)pab(O). ( 8 ) 

Since the QIP is monotonically nonincreasing under local 
CPTP maps acting on the system a, the function Q(p a b(t )) is 
monotonically nonincreasing with increasing time. Hence, 

^Q(Pab(t)) < 0 (9) 

dt 


-pit) = JXt)p(t), 

dt 


( 6 ) 


where J2(t) is the Liouvillian superoperator [3, 4], given by 

mp(t) = -iVH(t),p(t)] -^y,a)k(r)pa)A;w - 


-\A\{t)Ai{t),pif) }|. 


(7) 


holds for all t > 0 for a Markovian process. However, this may 
not be true for a non-Markovian process where the divisibility 
of the local CPTP map is violated. Therefore, j t Q{p a b(t)) > 0 
is a straightforward non-Markovianity witness. If we define 
D(t) — $- t Q(p a b(t)), then the non-Markovianity of the dynam¬ 
ical map can be quantified using the expression, 

A/q(Q) = max f D(t)dt, (10) 

Pad 0) JZ)(f)>0 
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where the maximization is performed over all sets of possi¬ 
ble initial system + ancilla states, p a b(0). The integration is 
extended over all time intervals for which D(1) > 0. Numer¬ 
ically, the final integration can be reduced to a summation of 
discrete sets of small interval integrals, 

A/a(Q) = max V f 'D k (t)dt. (11) 

Pab( '° ) k:D k (t )>0 ' 

The maximization over all possible initial states involved in 
quantifying non-Markovianity is certainly demanding. How¬ 
ever, starting with any chosen set of initial states, one can al¬ 
ways obtain lower bounds to the non-Markovianity measure, 
thus achieving a qualitative assessment of the non-Markovian 
character of the dynamics. 


III. QUALITATIVE ANALYSIS OF NON-MARKOVIANITY 

In order to qualitatively analyze the proposed non- 
Markovianity measure, A (q, and compare it with some of the 
other important measures in the literature, we now consider 
the paradigmatic model of a single-qubit dephasing channel. 
The Hamiltonian describing a single qubit interacting with a 
thermal reservoir is given by [1], 

'H = uj 0 a* + ^ cu/ata,- + + g*cr z d!), (12) 

i i 

where u>o is the qubit resonant transition frequency, a, (aj ) and 
at, are the annihilation (creation) operators and frequency of 
the /th reservoir mode, and g, is the reservoir-qubit coupling 
constant for each mode. The qubit dynamics resulting from 
the Hamiltonian in Eq. (12) is given by the differential equa¬ 
tion. 


pit) = y{t){crp{t)cr - p(t)), (13) 


where y(t) is the time-dependent decoherence rate, which can 
be determined from the spectral density (J(oj)) of the coupling 
constants [50, 51], and cr z is the third Pauli operator. The 
dynamical map for the dephasing channel (A(f)), on a single¬ 
qubit system (p a ) is given by 


Pa(t ) = A(f)p a (0) = 


P«Oo(0) Palo(0)r(0 
,PaOl(°) r (0 Pa 1 I (0) , 


(14) 


where T(f) = exp [-2 f Q ' y(t')dt'^, and p fl;y (0) are the elements 
of the initial system state, p fl (0). For a zero-temperature reser¬ 
voir with spectral density J{cS), the decoherence rate is given 
by the relation [1, 52], 


y(t) = 


f 




(15) 


while the ancilla is not subject to decoherence) optimized 
over all possible initial states p a b( 0). Such an optimization 
process is complicated and can be solved only for specific 
instances. Alternately, a lower bound on Nq can be ob¬ 
tained by considering the particular situation where the ini¬ 
tial composite system + ancilla state, p a b( 0), is maximally 
entangled, say p a bi 0) = |<D)(<1>|, where |<I>) is a Bell state, 

|<D> = (|00> + |11»/V2. 

The composite dynamical map is given by [28], 


Pabit) = (A a (0 <8> l*)|d>><0| 


1 

2 


i o o no 
oooo 
oooo 
T(0 0 0 1 


\ 




(16) 


and the non-optimized measure of non-Markovianity is given 
by the relation, 

Ni(Kit)) = f ow>() D(t)dt. (17) 

Atf(0)=l'F>CF| 

We need to calculate the QIP in the evolved state p a b{t) us¬ 
ing the measure <3, given by Eqs. (3)-(5). For the state p a b(t), 
given by Eq. (16), the maximum eigenvalue of the W matrix, 
A!" ax , is equal to 1 - T(f) 2 . Hence, 

Q(Pabit)) = vi - = no (is) 

D(t) = dQipabit)) = -2 T(0y(0- (19) 

Since T(0 > 0 for all f, it follows that D(t) > 0 when 
y{t) < 0, which is consistent with other well-established def¬ 
initions of non-Markovianity [6, 8, 28], The measure of non- 
Markovianity based on QIP for the initially maximally entan¬ 
gled, composite system + ancilla state, p a b, under single-qubit 
dephasing on the system a, is given by 

JVqWO) =-2 f Y(t)y(t)dt. (20) 

J y{t)< 0 

Interestingly, the quantification of non-Markovianity in terms 
of QIP (<3), for the paradigmatic single-qubit dephasing 
model and maximally entangled initial states, is numerically 
equivalent to the previously introduced measure in terms of 
the distinguishability of a pair of evolving states using the 
trace distance [8]. The distinguishability witness for non- 
Markovianity is closely associated with the backflow of in¬ 
formation from the environment to the system, which results 
in the increase of quantum correlations in the dephased bipar¬ 
tite state p a b as detectable through the QIP. Hence, the non- 
Markovianity defined in terms of the local quantum Fisher in¬ 
formation, whose minimization defines the QIP, exactly cap¬ 
tures the intrinsic backflow of information in the system- 
environment interaction. 

To further compare the measure Nq against other mea¬ 
sures of non-Markovianity, we consider specifically the above 
single-qubit dephasing model with an Ohmic reservoir spec¬ 
tral density, 



To analyze the non-Markovianity, we need to calculate the 
measure Nq, given by Eq. (10), for the composite system 
+ ancilla state, p a b (where the system undergoes dephasing 


J(a>) = aa>, 


( 21 ) 
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FIG. 2. (Color online). In this figure, we compare the measure Nq 
(= A( Q 0 ), with other quantifiers of non-Markovianity, for a single¬ 
qubit dephasing channel with an Ohmic spectral density, J(u>) = 
au> c ( a>/a> c ) s exp ( u>/u > c ). The blue dashed region, corresponding to 
y(t) > 0, is the Markovian regime. This corresponds to N = 0, for all 
known measures of non-Markovianity. The region y(t) < 0 denotes 
the non-Markovian regime and we plot the behavior of the measures 
based on QIP (Nq, green bar), quantum mutual information (TV/, red 
bar), and divisibility criteria ( N R , purple bar). We observe that the 
non-Markovian regime lies in the super-Ohmic region, 5 > 2. The 
behavior of the measures, Nq and Ni, is similar: Both indicators 
initially increase, followed by a decrease, with increasing S, with 
vanishing values for S >5. However, in contrast the measure Nr 
monotonically increases with S in the non-Markovian region. This 
implicitly shows that the measures based on QIP and mutual infor¬ 
mation are independent of the divisibility criteria. The inset figure 
shows the measures of non-Markovianity around the critical transi¬ 
tion parameter value, S = 2. We note that, for this model, the non- 
Markovianity measure based on distinguishability is identical to the 
one derived in this work. 


where S is the Ohmicity parameter, a is the dimensionless 
coupling constant, and a> c is the cutoff spectral frequency. For 
a zero reservoir temperature, relation ( 15) can be written as, 


7(f) = aa> c 


cos [5 tan '(We?)] ro(S) 
(1 + at 2 t 2 ) s/2 


( 22 ) 


where ToO) is the Euler Gamma function. It is known that 
the dephasing dynamical map corresponding to the spectral 
density J(a>) is divisible in the parameter range 0 < S <2 
[52]. Hence, the non-Markovian regime corresponds to the 
super-Ohmic parameter range S > 2, which can be experi¬ 
mentally obtained in ultracold systems utilizing control over 
atomic noise [18, 51]. To calculate Nq, we need to integrate 
over the relevant range of t, where y(t) < 0 and 'D(t) > 0. 

Figure 2 compares different measures of non-Markovianity 
based on QIP, quantum mutual information [28], distinguisha¬ 
bility [8], and divisibility criteria [6]. For a set of maximally 
entangled initial states, we observe that the measures based on 
the first three quantities behave in a similar manner as opposed 
to the non-Markovian measure based on divisibility. The mea¬ 
sures based on QIP and distinguishability are equivalent for 
this model. 


IV. FLOW OF QUANTUM INTERFEROMETRIC POWER 
IN NON-MARKOVIAN DYNAMICS 


We now consider the flow of the QIP (Q) in the dissipative dy¬ 
namics governed by the single-qubit amplitude damping chan¬ 
nel. The qubit dynamics can be modelled using the Hamilto¬ 
nian given by [1], 

V. = a>ocr 1 + 'Y_ i Uia]a t + ^(g;cr + a ; + g*cr~aJ), (23) 

i i 

where, cr + (cr ) is the raising (lowering) Pauli operator. The 
resulting dynamics is given by the differential equation, 

P = -*^[o-V~,p] + y(t) | o--pcr + - l - [crV~,p)j , (24) 


where p is time dependent. The functions s(t ) and y(t ) are de¬ 
fined in terms of the integro-differential equation for the time- 
dependent function, STt, given by 

j, = - f dTf(t-T)J„ (25) 

Jo 


where, s(t ) = -2Im^- and y(t) = -2Re^. The function ff t is 
characteristic of the nature of the environment used to model 
the local noise in the dynamics. It is defined in terms of the 
correlation function fit - r), that is derived from the Fourier 
transform of the spectral density of the environment, J(u>)\ 


fit - T) = J J(oj) (26) 

The dynamics of the system qubit ( p a ), is given by the dynam¬ 
ical map Pa(t) = Q.(t)p a ( 0), such that [1, 53] 

Pao 0 (0) + pan (0)( 1 - Jr) Pao 1 (0)JT, 


Pa(0 = 


(27) 


Paio(0)J7 Pan(°)J? > 


The single-qubit amplitude damping dynamics can be ex¬ 
tended to the system + ancilla bipartite system. In the com¬ 
putational basis, the evolved density matrix of the composite 
two-qubit system p a b(t) is given as follows (omitting the sub¬ 
script ab for simplicity). The diagonal matrix elements of p a b 
are given by [13, 54] 

Pn,ii(0 = Pii,ii(OX7mPio,io( 0 - Pio,io(OX7, 2 ; 

Poi,oi(0 = Poi,oi(0) +Pn,ii(0)(l _ jfy, 

Poo,oo(0 = 1 - (poi,oi(0+Pio,io(0+Pn,n(0)- 

The nondiagonal elements are given by 

Pii.toW =Pii,io(OX7?;pn,oi(0 =Pn,oi(OX7,; 

Pii.ooM =Pn,oo(OX7,;pio,oi(0 =Pio,oi(OX7,; 

Po 1 , 00 (f) = Poi,oo(0) + pn,io(0)(l _ Jj)\ 

Pio,oo(0 = Pio,oo(0)Jt, 

where p,- ; = p* ; ( i , j = 00,01,10,11). 
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Let us now consider a reservoir spectral density with a 
Lorentzian distribution [1,2], 


J(oj) = 


Tod 2 _ 

2tt[(w - wQ 2 + T 2 ] ’ 


(28) 


where lo c is the central frequency of the distribution and yo is 
the system-reservoir coupling constant. The spectral width of 
the distribution. A, is the inverse of the reservoir correlation 
time (t,. = 2). The system-reservoir coupling yo is related to 
the Markovian decay of the system, and is thus the inverse of 
the system relaxation time ( t s - -^). The Markovian nature of 
the dynamics is related to the strength of the system-reservoir 
coupling and the interplay of the system relaxation and reser¬ 
voir correlation times. For weak coupling, the relaxation time 
of the system is greater than the reservoir correlation time, 
t s > 2T r (yo < 4), and the dynamics is essentially Marko¬ 
vian. For, t s < 2 T r (yo > 4) or in the strong coupling regime, 
the dynamics is non-Markovian. Hence, the non-Markovian 
character of the considered dynamical map is ingrained in the 
behavior of the function 3t, which for a Lorentzian spectral 
distribution is of the form [1, 53], 


J, = e 


-(A-iS)t 



(d - iS) . , 

-sinh 



(29) 


where ij = yj(A- id) 2 - 2yod, and 5 (= (a>o - o> c )) is the 
system-reservoir frequency detuning. For the dynamics to 
be Markovian, in the weak coupling regime, the function J\ 
needs to have a monotonic decrease with time. For non- 
Markovian dynamics, the monotonicity of J~t does not hold, 
consistent with the breakdown of the divisibility of the dy¬ 
namical CPTP map. 

Now, let us consider our initial bipartite state to be a Bell- 
diagonal Werner state of the form. 


Pab( 0) = - 


I 4 + 2 ri(r ‘ 

V i=x,y,z 


(30) 


where the Werner parameter is given by, r x - —r y = r z -re 
[0,1]. The initial state p a b( 0) is maximally mixed for r - 0, 
and a maximally entangled Bell state for r = 1. The resulting 
time-evolved density matrix, under a single-qubit amplitude 
damping, can be written as 



f r + Jf 

0 

0 

2 rj, 

1 

0 

r-jf 

0 

0 

4 

0 

0 

r- + r + 3j 

0 


. 2 rJT 

0 

0 

r- + f J} 


(31) 


where, r ± = 1 ± r, and 3} = 1 -J~f. Using Eqs. (4), (5), 
and (29), the dynamical flow of QIP (Q) can be numerically 
evaluated for the density matrix given by Eq. (31). The non- 
Markovian character of the local dynamics of the bipartite 
qubit system ab can be analyzed by observing the evolution 
of Q(p n b(t)), with increasing time 1 of evolution. 

We begin by analyzing the case of a maximally entangled 
initial state p a b( 0), obtained by setting r = 1 in Eq. (30). If we 



FIG. 3. (Color online) The flow of QIP (<3) in a two-qubit system 
with maximal entanglement at f=0 under a single-qubit amplitude 
damping channel. The initial entanglement is set by setting r = 1 in 
the Werner state given by Eq. (30). Q(p a b ) is equal to \JJ ,|, as shown 
by Eq. (32). We observe the flow of Q in both Markovian and non- 
Markovian regimes of the dynamical evolution. The flow under the 
Markovian regime (red dashed line) corresponds to the ratio of the 
reservoir correlation to system relaxation, Ajy 0 = 10. Under the non- 
Markovian regime, the flow is shown for A/yo = 0.5 (blue dotted line) 
and A/yo = 0.1 (green solid line). The system-reservoir frequency de¬ 
tuning is set at, S = 0.01 y 0 . The backflow of quantum correlation, in 
terms of QIP, is observed by the nonmonotonic increase of Q during 
the evolution. 


consider/// to be a complex number of the form, 3t = a,+i/3 t , 
where a t , /?, 6 R V a,, /?, 6 [0,1], the maximum eigenvalue 
A™ ax of the matrix W, Eq. (5), is equal to (1 - aj -yS 2 ). Hence, 
the QIP Qipob) is given by, 

QiPabit)) = Vl - K ax = = \Jrl (32) 

Therefore, for a maximally entangled initial state undergoing 
a single-qubit amplitude damping, the nonmonotonic flow of 
the QIP in the non-Markovian regime is exactly governed by 
the nonmonotonicity of the function J~, (cf. [54]). Inter¬ 
estingly, \J,\ also measures the maximal trace distance be¬ 
tween a pair of system states and hence quantifies the non- 
Markovianity in terms of distinguishability for single-qubit 
amplitude damping channels [55]. We find therefore that the 
flow of the QIP Q(p a b(t )) is again closely related to the back- 
flow of quantum information in the non-Markovian regime. 

Figure 3 shows the flow of the QIP measure for an ini¬ 
tial maximally entangled state, given by Eq. (30), with r - 1. 
The Markovian and non-Markovian regimes of the dynamics 
can be studied in terms of the function J~ t , for a Lorentzian 
reservoir spectral distribution, as mentioned in Eq. (29). The 
Markovian regime corresponds to y (l /T < 0.5, as shown in the 
figure for A/yo — 10. The non-Markovian regime, correspond¬ 
ing to strong system-reservoir coupling yo /A > 0.5, is shown 
for A/yo = 0.1 and 0.5. The system-reservoir detuning is 6 - 
0.01 yo- The figure shows that the non-Markovian flow of QIP 
is nonmonotonic, with increase in Q during certain evolution 
times. The non-Markovianity can be numerically evaluated 
using the expression for A (q in Eqs. (10) and (11). 

We have seen that the flow of QIP for a maximally entan¬ 
gled initial state is determined by the nature of spectral dis¬ 
tribution and is equal to the integro-differential function \ J~ t \. 
However, the situation is not so straightforward if the initial 
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FIG. 4. (Color online). The flow of QIP (<2) for an initially mixed 
two-qubit system under a single qubit amplitude damping channel. 
The initial mixed state is obtained for r = 0.45, in the Werner state 
given by Eq. (30).The parameters pertaining to the Lorentzian reser¬ 
voir spectral distribution are set at A = O.Olyo, and S = O.OOlyo. 
The figure shows the non-Markovian evolution of Q(p a b ) (blue solid 
line), concurrence, C(p a t,) (black dotted line), quantum mutual infor¬ 
mation, Kfiab) (red broad-dashed line), and the scaled function |J7r|/2 
(brown dashed line). 


system + ancilla state is mixed. This can be obtained eas¬ 
ily from the Werner state in Eq. (30), by setting the Werner 
parameter r < 1. In such instances, the flow of QIP is still 
governed by the (non-)monotonic behavior of |J7?I- We ob¬ 
serve that the QIP (<3) has a dynamical behavior that quite 
closely replicates the flow of the function \SF t \, as compared to 
other measures such as the quantum mutual information [56] 
or the entanglement (quantified by the concurrence) [57]. Fig¬ 
ure 4, shows the dynamics of quantum information measures 
in the case of an initially mixed Werner state /»„/,. defined for r 
— 0.45, in the non-Markovian regime of the single-qubit am¬ 
plitude damping channel. The reservoir relaxation is set at 
A — O.Olyo, and the detuning is 6 - O.OOlyo. The figure shows 
specifically the flow of QIP, concurrence, and the mutual in¬ 
formation in comparison to the behavior of the scaled func¬ 
tion IJ7VI/2. For the considered model, the behavior of Q(p a b ) 
closely follows the nonmonotonic and discontinuous evolu¬ 
tion of the function J~ t . This is in contrast to entanglement 
and mutual information, which both evolve smoothly with 
time. Furthermore, entanglement decays quickly and vanishes 
for finite ranges of time (so-called entanglement sudden death 
[21]) and hence cannot qualitatively capture the backflow of 
quantum information in selected intervals of time. 

V. DISCUSSION AND CONCLUSION 

Non-Markovianity is an ubiquitous feature of quantum dy¬ 
namical maps, and is nowadays recognized as a resource for 
certain applications of quantum technology, such as metrol¬ 
ogy, cryptography, and communication [15-18]. The role 
of non-Markovianity in enhancing the robustness of quan¬ 
tum correlations in systems exposed to noisy environments 
has been studied by means of various quantitative approaches 
[6, 25, 28]. In this work, we have adopted the quantum 
interferometric power (QIP) [39, 40] as our reference fig¬ 
ure of merit to assess non-Markovianity of dynamical maps 


applied to a system coupled to an ancilla, which plays the 
role of a measuring apparatus for the operations occurring on 
the system. The QIP has been very recently acknowledged 
as a physically insightful, operationally motivated, and com¬ 
putable measure of quantum correlations of the most general 
kind, including and beyond entanglement [39, 58]. The QIP 
corresponds to the guaranteed metrological precision that a 
system + ancilla probe state enables the estimation of a phase 
shift on the ancilla part, in the black-box quantum metrol¬ 
ogy paradigm [39, 40]. When the system is subject to a 
non-Markovian evolution, the QIP between system and an¬ 
cilla (measured from the perspective of the ancilla) can un¬ 
dergo a nonmonotonic evolution, with revivals in time. We 
have shown that such nonmonotonic behavior is closely re¬ 
lated to, and can precisely capture, the backflow of informa¬ 
tion from environment to system which is a clear marker of 
non-Markovianity [8]. In operative terms, such a dynami¬ 
cal rise of quantum correlations translates into an increase 
of the guaranteed precision of phase estimation on the an¬ 
cilla, thanks to the non-Markovian noise affecting the system. 
While here we considered paradigmatic dynamical maps ap¬ 
plied to single-qubit systems only, it has been shown in [16] 
that non-Markovian noise affecting a register of n qubits can 
lead to an enhancement in the metrological scaling which 
is intermediate between the shot noise and the Heisenberg 
limit. It will be interesting to investigate how the measure of 
non-Markovianity proposed here in terms of QIP can be em¬ 
ployed to investigate the metrological scaling in the black-box 
paradigm for quantum metrology with multiqubit probes. 

In this work we have proposed to quantify non- 
Markovianity in terms of the nonmonotonicity of the QIP, 
similarly to previous proposals to quantify non-Markovianity 
in terms of the nonmonotonicity of entanglement or total cor¬ 
relation measures [6, 28]. By analyzing two simple mod¬ 
els of single-qubit noisy dynamics, we have shown how our 
measure reliably captures the non-Markovian regime, and is 
quantitatively more sensitive than measures based on entan¬ 
glement. Another advantage associated with the use of the 
QIP to characterize non-Markovianity, is that such a mea¬ 
sure also has been extended to Gaussian states of continuous 
variable systems, resulting in a computable and reliable mea¬ 
sure of quantum correlations in that relevant setting as well 
[43, 44]. In a subsequent work, it will be worth analyzing non- 
Markovianity in Gaussian dynamical maps [31, 32] in terms 
of the nonmonotonic flow of the Gaussian QIP [43]. The QIP 
has therefore the potential to offer a unified picture of non- 
Markovianity extending from qubits to infinite-dimensional 
systems. 

We hope that the present analysis can stimulate further re¬ 
search in order to pin down the relevance of non-Markovian 
dynamics in quantum information processing and in the de¬ 
scription and simulation of complex quantum systems in the 
biological, physical, and social domains [59]. 
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